Ground State Properties of a Tonks-Girardeau Gas in a Split Trap 
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We determine the exact many-body properties of a bosonic Tonks-Girardeau gas confined in a 
harmonic potential with a tunable (5-function barrier at the trap center. This is done by calculating 
the reduced single particle density matrix, the pair-distribution function, and momentum distribu- 
tion of the gas as a function of barrier strength and particle number. With increasing barrier height 
we find that the ground state occupation in a diagonal basis diverges from the \/N behavior that 
is expected for the case of a simple harmonic trap. In fact, the scaling of the occupation number 
depends on whether one has an even or odd number of particles. Since this quantity is a measure of 
the coherence of our sample we show how the odd-even effect manifests itself in both the momentum 
distribution of the Bose gas and interference fringe visibility during free temporal evolution. 

PACS numbers: 32.80.Pj, 05.30. Jp,03. 65. Ge 
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I. INTRODUCTION 

The last two decades have seen considerable experi- 
mental and theoretical activity and progress in the area 
of cooling and trapping of neutral atoms. Several crown- 
ing achievements have been produced, including the re- 
alization of Bose-Einstein condensation [ij , the creation 
of periodic arrays of single atoms [2] and the observation 
of superfluid phases in ultra-cold Fermi gases 01- While 
research in this area is interesting from a fundamental 
point of view, ultracold atoms are also well-suited can- 
didates to observe concepts and ideas in quantum infor- 
mation [3l . This is due to the fact that cold atomic sam- 
ples are often well isolated from the environment, while 
being highly controllable at the same time, which is of 
paramount importance if one wants to create and work 
with fragile many-body states. 

Two experimental advances have recently opened up 
the possibility to create and carry out experiments in 
strongly correlated quantum gases. The first one is due to 
optical lattices and atom chip traps, which can produce 
tightly confined potentials in selective directions of space. 
This allows one to limit the atom's degrees of freedom 
and thereby create effectively lower-dimensional systems 
Q. Secondly, by the using Feshbach resonances or 
by tuning of the effective mass of particles moving in a 
periodic potential Q , the inter-particle scattering length 
can be tuned to almost every value desired. 

Combining these techniques has permitted the experi- 
mental realization of atomicgases in the so-called Tonks- 
Girardeau (TG) regime [1, [J. A TG gas is defined to be 
a one-dimensional, strongly correlated gas consisting of 
bosons that interact via a hard-core potential 0, H, Q ■ 
In the limit of point-like inter-particle interactions Gi- 
rardeau found that such a model can be solved exactly 
by mapping it to an ideal, spinless fermionic system and 
he was the first to point out that a gas of strongly- 



interacting bosons can thereby acquire certain fermionic 
properties Q . More recently it was found that the above 
case is just a special case of a general mapping theorem 
between bosons and fermions in one-dimension, which 
can interact with finite strength 10]. 

In order to find the many-particle solutions of a given 
geometry for the TG gas (or for non-interacting fermions) 
one must first know the exact single particle eigenstates. 
However, since the list of exactly solvable single particle 
problems in quantum mechanics is limited, there is also 
only a small number of many-particle problems in the 
Tonks limit that can be exactly solved. Recently an ex- 
act solution for the experimentally important harmonic 
potential was found [ill . and here we will add another 
example to the list by describing a double well setting. 
As we aim for exact solvability, we have chosen the model 
of the (5-split trap [l^ , which comprises of a gas trapped 
in a harmonic oscillator potential split in the center by a 
point-like repulsive potential. In a recent paper the case 
of a boson pair was rigorously analysed in such a split 
trap for a range of interaction strengths up to and in- 
cluding the Tonks limit [lj|. It has also been suggested 
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in [15] that the (5-split trap could be used to excite dark 
soliton like structures in a Tonks-Girardeau gas. 

Investigating the model of the J-split potential can be 
justified in several ways. Firstly it can be seen as an ide- 
alized model of a realistic double well situation where the 
height of the barrier is related to the area of a physical 
potential. Comparison of our results with recent numer- 
ical simulations show that the qualitative behavior we 
find here persists for reasonable and realistic values of fi- 
nite sized splitting potentials [l^, [13, [3 • Alternatively, 
a point like potential can be a good approximation to 
describe a strongly localized impurity within the bosonic 
gas [l9l,f20i. 

Our main findings concern the coherence inherent in 
such samples. For the harmonically trapped case and in 
the thermodynamic limit it is known that the ground- 
state occupation scales as VN [21(]. However, for smaller 
samples it has been demonstrated that deviations from 
the ^/N behavior exist and that macroscopic coherence 
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effects can still be present [1^ . Therefore in this work we 
put particular emphasis on examining the ground-state 
occupation fraction as we increase the height of the cen- 
tral ^-barrier. 

The paper is organized as follows. In Sec. HI] we in- 
troduce the many-body Hamiltonian, describe the asso- 
ciated single particle eigenfunctions and eigenvalues and 
review the Fermi-Bose mapping theorem. Sec. IIIII ex- 
plores the physical many-body properties of our model 
by calculating the reduced single particle density matrix 
(RSPDM) and the pair-correlation function. In Sec. IIVI 
we investigate the influence of the splitting strength on 
the ground state occupation number and relate our re- 
sults for this behavior to two experimentally realizable 
quantities, namely the momentum distribution and inter- 
ference patterns. Finally, in Sec. |V]wc make concluding 
remarks. 



II. MODEL HAMILTONIAN AND THE 
FERMI-BOSE MAPPING THEOREM 

A. System Hamiltonian 

We consider a gas of N bosons trapped in a tight 
atomic waveguide. The waveguide restricts the dynamic 
of the gas strongly in the transversal directions, such that 
in the low-temperature limit we can restrict our model 
to the longitudinal direction only [23|. In this direction 
we then consider a i^-split harmonic potential such that 
at low linear density the many-particle Hamiltonian can 
be written as 



constant via 
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Here m is the mass of a single atom, uj the frequency of 
the harmonic potential and k is the strength of the point- 
like splitting potential, which is located at x = 0. Since 
we assume low densities, only elastic two-particle colli- 
sion have to be considered and we can restrict the inter- 
action potential to depend only on the relative distances. 
For bosonic systems at low temperatures the atomic in- 
teraction potential itself can be well approximated by a 
point-like potential 



V{\xi ~ Xj\) = giD5[\xi - Xj\) 



(2) 



where gm is the ID coupling constant. This approxi- 
mation is well justified for non-resonant situations and 
the only reminiscence of the exact potential is given by 
the three-dimensional s-wave scattering length, a^o- For 
positive values of a^u the interaction is repulsive and for 
negative values of 03^1 it is attractive. Finally, the scat- 
tering length is related to the one-dimensional coupling 
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where C is a constant of value C — 1.4603 . . . [231 • 



B. Eigenstates and eigenvalues of the (5-split trap 

The single particle eigenstates of the delta-split har- 
monic oscillator have recently been discussed in detail 
[2^ and we will briefly review them here for complete- 
ness. To do this we rescale the single particle part of the 
Hamiltonian ([1]) 



h 



2dx^ 



K,S{x) , 



(4) 



where all length are in units of the ground state size, qq — 
y^h/muj, and all energies in terms of the oscillator energy, 
huj. This leads to a new scaled length given by a; = xqq 
and scaled barrier strength given by k = (?jwao)~^K. For 
notational simplicity we shall drop the overbars on all 
scaled quantities and acknolwege that we are, henceforth, 
dealing in the scaled units just described. All units used 
in figure plots in this paper are also in terms of these 
scaled units. The time- independent Schrodinger equation 
for this system now reads 



hlpnix) = Enlpnix) 



(5) 



It is immediately clear that the odd eigenfunctions of the 
simple harmonic oscillator are still good eigenfunctions 
for the (5-split oscillator, as they vanish at the exact po- 
sition of the disturbance 



n — 1,3,5. 



(6) 



Here the Hn{x) are the n*"^ order Hermite polynomials 
and the Mn are the associated normalization constants. 
The corresponding energies are given by the eigenvalues 
of the odd parity states of the harmonic oscillator, £"„ — 
(n+i). 

The even eigenstates of the simple harmonic oscillator 
on the other hand have an extremum at a; = 0. They are 
therefore strongly influenced by the splitting potential 
and can be found to be ^] 
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where the U{a,b,z) are the Kummer functions The 
corresponding eigenenergies, i?„, are determined by the 
roots of the implicit relation. 
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Increasing the barrier height leads to an increase in the 
energy of the even eigenstates and in the limit of k = 
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oo each even eigenstate becomes energetically degenerate 
with the next higher lying odd eigenstate. 

Since we know the single particle eigenstates, we can 
build and solve the Slater determinant for a system of 
non-interacting fermions. Using the Fermi-Bose mapping 
theorem we can then calculate the many-particle bosonic 
wave- function from the fermionic result. 



C. The Fermi-Bose Mapping Theorem 

While the original Fermi-Bose mapping theorem onl 
related strongly interacting bosons to ideal fermions 
it was recently found that the mapping idea can be ap- 
plied to other systems as well [TO]- Here we concentrate 
on the situation relevant to our system, i.e. the Tonks 
limit of infinite, point-like repulsion {giD oo) between 
bosons. The main idea is that one can treat the inter- 
action term in eq. ([T]) by replacing it with the following 
boundary condition on the allowed bosonic wave- function 
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0, 



(9) 



for i j and 1 < i < j < N. As this is formally 
equivalent the the Pauli exclusion principle, one can solve 
for the associated ideal fermionic wave-function 
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det [V'„(xj)], 



(10) 



and calculate the bosonic solution from this by appropri- 
ate symmetrization 

*B = A{xi, . . . ,xn)^f{xi,X2, ■ . - ^xn), (11) 
where the unit antisymmetric function is given by Q 

(12) 



A= Y]. sgn(a;j - xj) 

l<i<j<N 



As we are only interested in the ground state, this last 
step simplifies to 



*s(a;i, xn) = l^pixi, xn)\ . 



III. MANY-BODY PROPERTIES 



(13) 



With the armory of the mapping theorem and the sin- 
gle particle states available, we are now in a position to 
calculate various ground state properties of the many- 
particle state as a function of increasing particle number 
as well as varying the strength of the central S barrier. 



A. Reduced single particle density matrices 

As laid out in the previous section, the Fermi-Bose 
mapping theorem allows for the calculation of the exact 



many-particle wave-function. For the case of an infinitely 
strong barrier this was recently done analytically in [l3| 



*s = 



C 



N/2 2 {N/2M/2) 

i J 0",fe)=(ij+i) 

(14) 

where C is the normalisation constant. Although this 
function fully characterizes the state of the system, other 
quantities can be more useful for obtaining characteristic 
properties of many-particle systems. We therefore pro- 
ceed to calculate the reduced single particle density ma- 
trix (RSPDM), from which expectation values of many 
important one-body physical observables such as the mo- 
mentum distribution or the von Neumann entropy can 
easily be obtained. 

For the bosonic gas, the RSPDM is defined as 



p{x,x')^ / 'if*B{x,X2, ■ ■ ■ ,Xn)^ 
J —oo 

b{x' ,X2, ■ ■ ■ ^xn) dx2 ■ . -dxN , (15) 

and we choose its normalisation to be given by 
J p{x,x)dx = N. While for the simple harmonic oscil- 
lator this integral was solved analytically by Lapeyre et 
al. (25j . we have to resort to a numerical evaluation for 
finite values of n. Here we use an algorithm recently pre- 
sented by Pezer and Buljan [26|| that allows for effective 
calculation of eq. (fTS]) for large numbers of particles. In 
fact, our particle number is only limited by the numer- 
ical instabilities when calculating higher order Kummer 
functions. 

The RSPDM expresses self correlation and one can 
view p(x,x') as the probability that, having detected 
the particle at position second measurement, im- 
mediately following the first, will find the particle at the 
point x'. Classically, p{x,x') = 5{x ~ x'). The RSPDM 
for a twenty particle Tonks gas is shown in Fig. [T] In 
the unsplit case (k = 0) we see that most of the density 
is concentrated along the diagonal. Increasing k intro- 
duces a gap around the position of the splitting potential 
a; = x' = 0, along with a reduction of density in the 
off diagonal regions. This is due to the suppression of 
tunnelling from one side of the barrier to the other as 
the systems eigenstates become doubly degenerate in the 
K — > cx) limit. 

The single particle density can be calculated from the 
RSPDM by taking x = x' 



p{x) = N 



+ CXD 



|*s(a:,a;2 . . .,XN)\'^dx2 



. dx 



N 



(16) 



and is therefore simply given by the diagonal of p{x, x'). 
In Fig. [2]we show this quantity for several different values 
if K. In contrast to the self correlations, one can see 
that the splitting effects the density only very locally. 
The results found in this case completely agree with the 
results presented in ^l3| , where they were calculated in a 
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FIG. 1: (Color online) Reduced single particle density matri- 
ces p{x,x') for a Tonks gas of twenty particles for splitting 
strengths n — (left) and k = 20 (right). 



3 1 1 3 1 1 3 




-10 10 -10 10 -10 10 



XXX 

FIG. 2; (Color online) Single particle density for a Tonks gas 
of twenty particles for increasing barrier height. 



different way by summing up the single particle densities 
and symmetrising. 

To gain further understanding we next calculate the 
eigenvalues and eigenfunctions of the RSPDMs given by 

/oo 
p{x, x')<j)j{x')dx' — Xj(j)j{x) . (17) 
'OO 

The eigenfunctions (x) are known in theoretical chem- 
istry as natural orbitals and their associated eigenvalues, 
\j , represent the occupation number of each orbital. The 
first three lowest energy natural orbitals for a twenty par- 
ticle gas with a splitting potential of height k = 20 are 
displayed in Fig. [S] One can see the point-like distur- 
bance introduced by the (5-function into the symmetric 
orbitals 0o and 02, while the anti-symmetric (f)i is un- 
changed from the k = case (natural orbitals for the 
K = case are displayed in [T^). 

These states and their occupation numbers will be used 
in section HVl to calculate the reciprocal momentum dis- 
tributions and the ground state occupations. 



B. Pair distribution functions 

The pair-distribution function D{xi,X2), is a two- 
particle correlation function that describes the probabil- 
ity to measure two atoms at two given positions at the 
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FIG. 3: (Color online) The first three natural orbitals for 
a twenty particle Tonks gas in a 5-split trap with splitting 
strength k = 20. 

same time. It is defined in the following way 

/ + 00 
\'i'B{xi,X2 . ■ .,XN)\'^dX3 . . .dxN 
-OO 

(18) 

N-1 
0<n<n'<N-l 

(19) 

where Afjj = N{N — 1). Since the terms with n = n' 
in eq. (llOp vanish, and we can rewrite it in the following 
form 

D{xi,X2) = p{xi)p{x2) - |A(a;i,a;2)P , (20) 

which is dependent only on the single particle density 
and the correlation function |A(a;i,a;2)| defined by 

|A(X1,X2)|= rn{xi)i^n{x2) ■ (21) 

0<ri<ri'<Ar-l 

The pair distribution functions for samples consisting of 
N = 5, 10 and 30 particles and for barrier heights of 
K = 0, 1 and 10 are shown in Fig. [4l The first striking 
feature inherent to all situations is the absence of any 
probability along the diagonal, which is due to the im- 
penetrable nature of the atoms. As the splitting strength 
is increased one observes the absence of probability for 
a joint measurement along the cross defined by the lines 
xi = and a;2 = 0, which is again due to the central 
position of the i5-splitting. 

IV. GROUND STATE OCCUPATION 
NUMBERS AND COHERENCE EFFECTS 

The fraction of particles that are in the (I)q{x) orbital 
/ is related to the largest eigenvalue Aq of the RSPDM 
by / = Therefore, in analogy to the macroscopic 
occupation of a single eigenstate in a Bose-Einstein con- 
densate, this orbital is sometimes referred to as the BEC 
state and the quantity Aq hence acts as a measure of the 
coherence in the system. Recently Forrester et al. have 
shown that, as one increases the particle number, Xq/N 
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FIG. 4: (Color online) Pair distribution functions D{xi,X2) 
for a Tonks gas of A*' = 5, 10 and 30 particles in 5-split trap for 
splitting strengths k = 0, 1 and 10. In each plot the horizontal 
and vertical axes run from —10 to +10 in scaled units. 



A. Occupation numbers 

The fractional ground state occupation Xq/N is dis- 
played in Fig. [5] as a function of particle number for 
different heights of the splitting potential. The dotted 
line corresponds to the unsplit trap (k = 0) and agrees 
with previously published results [2l|, [23| • When increas- 
ing the magnitude of the splitting one notes a strikingly 
different behavior for the values of Xq/N for odd and 
even particle numbers. For even particle numbers we 
find the coherence decreased as compared to the k = 
case, whereas it is essentially unchanged when the parti- 
cle number is odd. The effect becomes more pronounced 
as the splitting strength is increased and it damps out 
as the particle number is increased. This behaviour is 
directly related to the magnitude of the RSPDM in the 
off-diagonal quadrants, sgn(a;) ^ sgn(a;') and therefore a 
clear signature of the coherence inherent in the system. 
To see this behaviour manifest itself in an experimentally 
realizable quantity we will next calculate the momentum 
distribution as well as the visibility of interference fringes 
in the free evolution of the gas when released from the 
trap. 



B. Momentum distributions 
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FIG. 5: Plot of the fractional ground state occupation f ~ ^ 
versus particle number for a TG gas in a 5-split trap with 
increasing splitting strength. 



tends toward l/y/ N in the harmonically trapped case 
[2l| . Here we will show how the introduction of a central 
barrier affects this 1 / \/N behavior in a dramatic way. 



For a harmonic trap the relationship between the mo- 
mentum distribution and coherence properties of the TG 
was recently studied by Minguzzi and Gangardt psl ] . The 
momentum distribution, n(fc), can be calculated from the 
reduced single-particle density matrix 



n{k) EE (2?:)- 



+ OC /' + OC 



p(x,a;')e-''=("'-^')da;dx' , 



(22) 

and is normalised as n{k)dk ~ N . Equivalently it 
can be obtain by considering the eigenstates of p{x,x'). 
Using a discretised form for the quadrature then allows 
one to rewrite the integral equation as a linear algebraic 
equation 



i(fc) = 5]A.|M.(fc)|' 



(23) 



where ^i{k) denotes the Fourier transform of the natural 
orbital 4>i{x), 



A*i(fc) = 



1 



27r J- 



+ 00 



(24) 



The central peaks of the reciprocal momentum distri- 
bution for N — 5 and = 6 particle samples are shown 
in Fig. [5] for different values of the splitting strengths. 
Fig. [7] shows the same quantity for gases with A'' = 19 
and N — 20 particles. While the states with odd parti- 
cle numbers are clearly less affected by the barrier than 
the states with even particle number, one can see the 
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FIG. 6: (Color online) Central peaks of the momentum dis- 
tributions n{k) for TG gases consisting of A'^ = 5 (left) 
and A'^ = 6 (right) particles, for values of splitting strength 
K = 0, 10 and 20. 
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FIG. 7: (Color online) Central peak of the momentum distri- 
bution n{k) for TG gases consisting of A'' = 19 (left) and 
N = 20 (right) particles, for values of splitting strength 
K = 0, 10 and 20. 



emergence of bi-modality at the neck of the peaks in the 
plots on the left hand side. It stems from the interfer- 
ence of particles on both sides of the splitting potential. 
For even particle number the introduction of the central 
splitting significantly broadens the momentum distribu- 
tion lowers it's peak. This is in agreement with analyti- 
cal results we have found earlier for the special case of a 
two-particle Tonks molecule [3] and indicates a loss of 
coherence within the sample. The form of the momen- 
tum distributions displayed in Fig. [HI and Fig. [7] is deter- 
mined by Eq. 1231 In this equation Aj acts as a weight 
on the contribution of the Fourier transform of an in- 
dividual natural orbital. It is interesting to see which 
natural orbitals are responsible for the broadening of the 
distribution in the even particle case. In Fig. [5] the first 
six Aj's are plotted as a function of k for iV = 19 and 
TV = 20 particles. The difference between the odd and 
even samples can be clearly seen. In the TV = 20 case 
we see a degeneracy of consecutive eigenvalues occuring, 
this is due to the symmetry of having ten particles on the 
right of the barrier and ten on the left. The degeneracy 
explains why we have a noticable change in the respec- 
tive momentum distribution. In the = 19 case no such 
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FIG. 8: Magnitude of the first six eigenvalues \j of p{x, x') 
for A'^ = 19 and A'^ = 20 particles as a function of the splitting 
strength k. The lines of greatest magnitude represent Ao and 
decrease in the order Ao, Ai, A2, A3, A4, A5. 



degeneracy occurs due to fact that one of the particles 
is spatially delocalised over both traps, this also explains 
the emergence of bi-modality in the neck of the odd dis- 
tributions. This behavior was found to be consistent for 
all particle numbers studied, up to 30 particles. 



C. Interference fringes 

The Fermi-Bose mapping theorem holds equally well 
for time dependent ground-state wave- functions and we 
will study the time evolution of the many-body quan- 
tum state after removal of the external potential. This 
situations is similar to one studied by Girardeau and 
Wright |11|, who considered splitting and recombining 
a gas within an external trap. In our system splitting 
the sample is inherent in the Hamiltonian and we will 
pay special attention to effects stemming from different 
particle numbers. Starting off with the gas confined in 
the 5-split trap with large splitting amplitude, we look 
at the time evolution of the single particle density p(x, t) 
as both, the trap and the central splitting, are turned 
off and the gas undergoes free temporal evolution. Dur- 
ing this both halves of the trap start overlapping and the 
densities for samples with N = 5 and N = 6 particles and 
with iV = 19 and iV = 20 particles are shown in Figs. [9] 
andfTPl respectively. In both cases one can clearly see a 
more distinct interference pattern emerging for odd par- 
ticle than for even particle samples, indicating that odd 
states carry larger coherence. In agreement with Fig. [S] 
the coherence effect is less pronounced at larger parti- 
cle numbers but still clearly visible as can be seen from 
Fig.M 

It is well worth to point out that as Figs. [9] and [10] show 
only density distributions, the results apply equally to a 
gas of spin polarized fermions. In this case, however, 
the interpretation of the variations of coherence with re- 
spect to particle number is straightforward, as the split- 
ting introduces non-smooth kinks in the single-particle 
wave- function which create the Slater determinant. 
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FIG. 9: (Color Online) Free temporal evolution of the single 
particle density, p{x,t) for a A*" = 5 and A'' = 6 particle TG 
gas initially confined in a 5-split trap with k = 100. 
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controlled way due to its analytic accessibility. Identi- 
fying and describing new potentials that take advantage 
of this is therefore of large importance. In the present 
work we have undertaken a thorough investigation of the 
many-body properties of the Tonks-Girardeau gas in a 
(5-split trap. We have calculated the RSPDM as well as 
the pair distribution function for various different split- 
ting magnitudes and particle numbers and identified the 
basic physical behaviour shown by the system. 

From the RSPDMs we were able to study coherence 
properties of the gas by determining the behaviour of 
the ground state eigenvalue, Aq, as a function of parti- 
cle number. Our results show that odd and even particle 
number samples obey different scaling laws, with the odd 
number samples remaining more coherent or less sensi- 
tive to the central splitting. The effect becomes less pro- 
nounced as one goes to larger particle numbers. To show 
how this effect manifest itself in different and experimen- 
tally observable quantities, we have studied the momen- 
tum distribution and interference experiments. For the 
momentum distributions we found that for odd particle 
numbers the sharp peak around momentum zero is rela- 
tively insensitive to the different strength of the splitting. 
For even particle numbers, however, the distributions are 
lowered and widened with increasing k, demonstrating a 
loss in coherence. The simulations of the interference 
experiments for odd and even samples showed a larger 
visibility occurring for odd particle samples. This is in 
agreement with the other quantities that the odd number 
samples are more phase coherent. 



FIG. 10: (Color Online) Free temporal evolution of the single 
particle density p{x,t) for a = 19 and A'' = 20 particle TG 
gas initially confined in a 5-split trap with n = 100. 



V. CONCLUSIONS 

Tonks gas has over the previous years shown to be 
an exciting and rich system to study new physics in a 
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